Vortex relaxation and coupling in superconducting heterostructures studied by STM 
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In a sandwich consisting of two superconducting films, one weakly pinning and one strongly 
pinning, the vortex positions in both films are determined by the strongly pinning material and the 
vortex lattice is disordered in both films. We used (strongly pinning) NbN and (weakly pinning) 
a-Mo2.7Ge and studied, by directly imaging the vortex core positions with a scanning tunnelling 
microscope, how this disorder is restored with increasing thickness of a-Mo2.7Ge layer or when the 
interface is made insulating. For clean interfaces we find that the first reordering of the vortex 
lattice is found at a layer thickness wich is compatible with the first bending mode of the vortex 
lines. Making the interface insulating we find that order is restored quickly. We argue that this is can 
be understood from the competition between the Josephson force working on the vortex segments 
on the one hand, and the elastic restoring forces inside the weakly pinning layer on the other hand. 

PACS numbers: 74.25.Qt,74.50.+r 



I. INTRODUCTION 



In the investigation of the properties of vortices and 
vortex lattices in type-II superconductors, direct imag- 
ing of the vortex positions at the surface of the material 
can be a useful tool. For magnetic fields of practical in- 
terest the distance between vortex lines becomes small 
(roughly 50 nm at 1 T) and Scanning Tunnelling Mi- 
croscopy (STM) gains an advantage over magnetic imag- 
ing techniques since STM is only sensitive to the vor- 
tex cores. A disadvantage of the technique is that the 
sample surface needs to be clean and conducting. For 
crystals, this can be done by breaking or cleaving, which 
has been successful for such materials as NbSe2 [ll, Hj] or 
Bi2Ca2SrCu208-(5 0,0]. For thin films no such simple 
preparation method exists and since surfaces easily oxi- 
dize, reports on vortex imaging on superconducting films 
are relatively few. One way to overcome this problem 
is to protect the film surface with a thin capping layer 
of a noble metal such as Au. As we showed before, this 
works well on films of a-Mo2.7Ge, a weakly pinning su- 
perconductor, capped with 3 nm Au and clear images 
of the vortex structure could be obtained by mapping 
the proximity-induced gap at the Au surface [5]. The 
Au layer should be flat, which means that Au should 
be wetting the surface during deposition, and the inter- 
face should be transparent. Neither condition is always 
met. For instance, although a reasonably flat Au film 
forms on films of strongly pinning NbN, spectroscopy 
shows very unclear vortex images, presumably because 
the rough NbN surface causes a bad NbN / Au interface. 
This can be overcome by depositing a thin (20 nm) film 
of a-Mo2.7Ge followed by the Au protection layer. Clear 
vortex contrast could again be obtained, and showed the 
NbN vortex lattice (VL) to be completely disordered [5[. 
This use of a a-Mo2.7Ge/Au template layer assumes that 
the vortex positions as fixed by the NbN layer are trans- 
ferred in true fashion to the Au surface, which should be 
the case when the a-Mo2.7Ge layer is too thin to yield ap- 



preciable bending or rearrangement of the vortex cores. 
In this work we further address this issue, which consists 
of two parts. The first involves the vortex- vortex interac- 
tions within the template layer, which might lead to reori- 
entation of the vortices in that layer or, more specifically, 
to the recurrence of the well-ordered a-Mo2.7Ge VL. The 
second question is with respect to the coupling of the vor- 
tices in the two superconducting layers. If the interface 
has a metallic character, the proximity effect will play a 
major role in coupling the vortices; if the surface of the 
host superconductor is insulating, as is often the case in 
oxide materials, the coupling may become of the Joseph- 
son type, or merely electromagnetic, and this can allow 
freedom for the vortex in the template to move with re- 
spect to the one in the host. 

The structure of the paper is as follows. We first in- 
troduce the experimental techniques and the data anal- 
ysis process, in particular the use of the autocorrelation 
function. Then we address the question of the length 
scale on which the disorder injected by the NbN layer 
is healed in the a-Mo2.7Ge template layer, by measur- 
ing samples with different (template) thicknesses dt of 
the a-Mo2.7Ge-layer. We show data on the variation of 
the autocorrelation functions as function of dt^ and we 
discuss them by considering the bending modes of vor- 
tex lines in a-Mo2.7Ge. Then we present experimental 
results on the issue of coupling, by studying heterostruc- 
ture samples with a badly conducting or insulating inter- 
face layer. In the discussion, we consider the strength of 
the electromagnetic and Josephson coupling forces which 
couple the vortex segments in the two superconducting 
layers. This coupling is in competition with the inter- 
vortex forces within the template layer, which tend to 
reorient the vortices to their triangular arrangement. 
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II. EXPERIMENTAL DETAILS AND DATA 
PROCESSING 



In this work a-Mo2.7Ge is used as a template layer be- 
cause of its suitable properties. With a Ginzburg-Landau 
parameter k ^ 100 it is a strong type-II superconduc- 
tor, while its amorphous structure makes it a relatively 
weakly pinning material. Furthermore, its large Lon- 
don penetration depth Xl yields a small tilt modulus 
C44, which means that the VL retains a two-dimensional 
(2D) character up to quite large thicknesses, as will be 
discussed below. Finally, amorphous materials show very 
flat surface characteristics. The a-Mo2.7Ge layers were 
rf-sputtered from a composite Mo/Ge target in Ar at- 
mosphere. The film thickness ranged from 25 nm to 
100 nm. When necessary, we denote the thickness by a- 
Mo2.7Ge(x), with x the thickness in nm. During sputter- 
ing the substrate was water-cooled and after each minute 
of sputtering, and interruption of 1 minute was made to 
prevent heating and possible recrystallization of the films. 
The sputtering rate was chosen as ~ 8 nm/minute. The 
superconducting transition temperature Tc was 5.7 K for 
the 25 nm films, increasing up to 6.3 K for dt = 50 nm, 
after which Tc became constant. The critical current for 
onset of vortex motion for a 500 nm thick film was found 
to be or the order of 10^ A/m^. Also relevant is the value 
of the critical field Bc2 at the measurement temperature 
of 4.2 K, which was about 4.2 T for dt = 25 nm, in- 
creasing to 5.2 T for = 50 nm and beyond. The STM 
measurements were performed in fields up to 1 T, which 
means maximum reduced fields b = B/Bc2 ^ 0.2 - 0.3. 
The strongly pinning NbN films (Jc > 10^ A/m^) were 
reactively rf-sputtered from a Nb target in an Ar/N2 at- 
mosphere on Si substrates. Throughout the study, we 
used a thickness of 67 nm. The films are polycrystalline, 
as is known from high resolution transmission electron 
microscopy studies [6|. The rms-roughness of the films 
was determined by AFM to be 1.5 nm, measured on a 
fresh NbN surface area of 1 x 1 /im. The films showed 
a Tc of 11.4 K and a superconducting coherence length 
^ = 4.7 nm. 

Protective Au films were rf-sputtered in two steps. First, 
about 2 nm was sputtered in a pure Ar atmosphere; next, 
oxygen (partial pressure 5 %) was mixed in the sputter- 
ing gas, and another layer of 2 nm was sputtered. This 
02-addition promotes the formation of smaller Au grains 
which improves the resolution of the measurements. 

The STM-setup was the same as in a previous report 
[2], with the sample immersed in liquid He after it was 
mounted in air. In this procedure field-cooled measure- 
ments are not possible, since the level of the liquid is 
higher than the magnet. All data were therefore taken 
by zero-field-cooling the sample, and then increasing the 
magnetic field. For all experiments we employed full 
current (I)- voltage (V) spectroscopy, with a single curve 
consisting of 200-300 points. The voltage range was cho- 
sen such that a good fit could be made to the Ohmic resis- 
tance outside the gap region, in practise typically ±4 mV. 
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FIG. 1: Averaged normalized derivatives dl/dY for two sam- 
ples, acquired at T = 4.2 K and B = 0.3 T. (a) Sample 
NbN/MoGe(50)/Au(6); (b) Sample NbN/MoGe(50)/Au(10) 



Several curves were taken per pixel, and images typically 
contained 128 x 128 pixels. The total acquisition time 
for such an image was of the order of 30 minutes to 1 
hour, depending on the surface quality. The maximum 
scan size, for the given grid, is determined by the vortex 
core size. The vortex signature is visible in a region with 
a diameter of roughly 4 ^, and needs about 3 pixels in this 
region for a proper determination. This yields 7 nm per 
pixel, and a scan size of maximally 0.9 /im. The vortex 
positions were found by taking the ratio of the zero-bias 
and high-bias slopes, as obtained by numerical differenti- 
ation. Theoretically, this is a value between (full gap) 
and 1 (no gap). In practise the gap is not well devel- 
oped due to the finite temperature and the proximizing 
Au layer and in the gap regime the ratio is of the or- 
der 0.8. The resulting range (0.8 to 1) is then translated 
into a grey scale. Fig. [T] shows the effect of the thickness 
of the Au layer by comparing the gaps for two samples 
NbN/MoGe(50)/Au, where the first one has a Au layer 
thickness of 6 nm, while the Au layer of the second one 
is almost two times thicker. 

A useful way of determining the amount of correlation 
of vortex positions is to compute the two-dimensional 
autocorrelation function G{v) for a given image, defined 
as G(r) =< p(ri)p(ri + r) >. The function computes 
the overlap of the image with a copy of itself which is 
displaced by a vector r. In the actual calculation, the 
result is normalised in order to correct for the decreasing 
amount of datapoints available with increasing r. Fig. [2] 
gives the result of this procedure for a VL imaged on a a- 
Mo2.7Ge-film of 48 nm thick at an applied field B of 0.8 T. 
Triangulation (Fig. [2)3) shows that the VL is quite free of 
defects. The autocorrelation map of G(r) (Fig.lSb) shows 
very little blurring towards the edges, and makes evident 
that the VL has both high translational and rotational 
bond order. By making a cross-section along a direction 
from the center of the image through one of the points 
of the hexagon, the loss of correlations can be quantified. 
(Fig.[2]i) shows such a cross-section, plotted in units r/ao 
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FIG. 2: (a) Vortex lattice (VL) image on a sample 
MoGe(48)/Au, acquired at T = 4.2 K and B = 0.8 T. The 
scan size is 1.2 /jui x 1.2 /xm. The raw data were smoothed 
with a Gaussian filter of 1 pixel width; (b) Triangulation rep- 
resentation of the VL. The lines represent bonds, the symbols 
denote deviations from the 6-fold coordination of the perfect 
trangular lattice (triangles (squares) : 7-fold (5-fold) coor- 
dination); (c) Normalized autocorrelation map. (d) Cross- 
section of the autocorrelation map along the white line in (c) . 
The inset gives the central part of the 2D-Fourier transform 
of the direct-space image given in (a) 

and normalized to the value at r/ao = 0. Here, ao the in- 
tervortex distance given by ao = 1.07\/^/B. It appears 
that Tc (somewhat arbitrarily defined as the displacement 
vector where the autocorrelation function does not show 
lattice periodicity anymore) is larger than 16 ao 

In contrast, the VL in NbN is almost completely 
disordered, as shown in Fig. [3] for a 50 nm film with 
a Au capping layer. In this case, we find Tc ~ 1.5ao 
and almost no orientational correlation. Note also the 
relatively poor quality of the vortex image, due to 
the fact that no a-Mo2.7Ge-template layer is present. 
Next we shall investigate how this disorder propagates 
through an intervening a-Mo2.7Ge layer. 



III. HEALING OF THE VORTEX LATTICE 

A. Experimental data 

In order to investigate on what length scale the VL dis- 
order starts to approach the a-Mo2.7Ge hexagonal lattice 
again, we recorded vortex images at three field values 
(0.2 T or 0.3 T, 0.7 T or 0.8 T, and 1 T) for NbN/a- 




B = 0.2T scansize = 419 nm B = 0.5T scansize = 419 nm 




B = 0.7T scansize = 419 nm B = 0.7T scansize = 419 nm 




FIG. 3: Compilation of results of vortex lattice (VL) imaging 
of a sample NbN(50)/Au. The middle panel in the left-hand 
corner shows a VL image acquired at T = 4.2 K and B = 
0.7 T. The adjacent panel shows the normalized autocorre- 
lation map computed from theses data. Other panels show 
autocorrelation map at fields of 0.2 T, 0.5 T and 0.8 T. The 
lower right-hand panel shows cross-sections of the autocorre- 
lation images taken at 0.2 T and 0.8 T. 



Mo2.7Ge/Au samples with dt = 0, 25, 50 and 100 nm. 
In Fig. m we show the autocorrelation (acr) results com- 
puted from the raw vortex images as measured in either 
B = 0.7 or B = 0.8 T, for the samples with dt = 0, 25, 
100 nm, together with one real image for dt = 50 nm as 
a sample of the data quality. Since the contrast and the 
noise-level of the data are affected by the varying exper- 
imental conditions the relevance lies in the the period- 
icity of the structures. The gray-scale of the acr-images 
is chosen such that the structure of the data is clearly 
visible. From the data it appears that for the thinnest 
a-Mo2.7Ge layer, the vortex positions are basically un- 
changed. The first ring at r = ao stands out more clearly, 
but that is simply due to the better data quality. There 
is also some increased intensity at the hexagon corners, 
indicating the onset of orientational order. On increasing 
the thickness of the a-Mo2.7Ge layer we find that both 
the transversal and orientational order develops more and 
more. This is better followed in the cross-sections of the 
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MoGe(25) MoGe(50) 
B = 0.8 T - scan size = 627 nm B = 0.7 T - scan size = 471 nm 




IVIoGe(IOO) IVIoGe(IOO) 
B = 0.8T - scan size = 374 nm B = 0.8 T - scan size = 374 nm 

FIG. 4: Compilation of autocorrelation functions resulting 
from vortex lattice (VL) imaging of samples NbN(67)/a- 
Mo2.7Ge(dt)/Au, with dt = 25 nm, 50 nm, 100 nm acquired at 
T = 4.2 K and fields of either B = 0.7 T or 0.8 T as indicated. 
The lower right-hand corner shows real space data. 
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FIG. 5: Cross-sections from the autocorrelation maps 
of vortex lattices measured on samples NbN(67)/a- 
Mo2.7Ge(dt)/Au, with dt = nm, 25 nm, 50 nm, and 100 nm 
acquired at T = 4.2 K and fields of either B = 0.7 T or 0.8 T 
as indicated. 



acr-images presented in Fig. [5l Clearly, the reordering 
occurs but slowly. For dt = 25 nm, Tc/ao ^ 2 at all 
fields. For 100 nm, this has increased to about 4 (in the 
0.8 T data) or maybe 6 (in the 1 T data). Although the 
signature for healing therefore is unequivocally present, 
the order at a a-Mo2.7Ge thickness of 100 nm is still far 
less than in the unconstrained vortex lattice. 



B. Bending modes in a-Mo2.7Ge. 

In order to understand this length scale, we next con- 
sider the length over which a vortex line can bend. This 
involves both tilting the line and shearing it with respect 
to the whole VL. Starting from a 2D collectively pinned 
VL, we can estimate the film thickness where the first 
bending starts to be relevant. The extent to which de- 
formations can develop in the transverse and longitudinal 
direction depends on the balance of the tilt and shear en- 
ergy densities Eg^t which are given in the elastic limit by 
Eg = ^cqq{u/s)'^ for shear deformations (with cqq the 
shear modulus of the VL) and Et = |c44(ii//)^ (with 
C44 the tilt modulus of the VL) for tilt deformations of 
a vortex in a medium where s and / respectively are the 
transversal and longitudinal axis of the cigar-shaped sur- 
face with constant displacement u (with <C ao in order 
to be in the linear regime). The energy balance leads to 
a relation for the ratio between s and /: 



/C44(fe±,fe||) 
C66 



(1) 



This means that if we set the typical length scale for lon- 
gitudinal deformations to the film thickness (/ = d), the 
typical transverse length scale s is set by the square root 
of the ratio of the elastic constants. The smallest trans- 
verse length scale is simply given by 5 = ao- The thick- 
ness dmax where these deformations start to be relevant 
can be obtained by neglecting the very small dispersion 
of Cqq and taking for C44 = C44{k±^ ku): 



C44(0) 



C44(0) 



X'ikl^k^) 



x^ki^k^y 



(2) 



where 044(0) = B'^/jj^q; in the last step we made use of 
the fact that A ao, as well as that b < 0.3 (making a 
field correction for A unnecessary). Since A is very large, 
we use the non-local expression. The wavevector k± for 
tilt deformations runs from to Kq^ the Brillouin zone 
boundary. In a circular approximation the Brillouin zone 
boundary is given by Kq = (47r^5/<l>o)^^^. In principle 
one should integrate over all the values for k±^ but from 
eq. [2] we find that the large /c-contributions dominate, 
leading to the approximation: 



C44{k_i,k\\) 



C44(0) 



(3) 



In case of a given vortex length d the wavevector /cy for 
tilt deformations is quantized: k\\ = nn/d with n G N. 
Since we are looking for the case that the vortex line 
just starts to bend, we are interested in the case where 
the n = 1 mode becomes relevant, which leads to the 
expression 



C44{k±, k\\) 



C44(0) 



(4) 
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FIG. 6: Dependence of dmax, the maximum film thickness for 
which (tilt) deformations are absent in a-Mo2.7Ge-films, on 
the reduced magnetic field b 



It can easily be seen that the relevance of higher order 
terms goes down with n~^. Together with Eq. [T]we have 
an expression for the maximum film thickness dmax foi" 
which deformations of the vortex line by bending are ab- 
sent. The dependence of dmax on b is plotted in Fig. [6] 
using the parameters for a-Mo2.7Ge. We find that in 
the experimental regime of 0.01 < b < 0.3 for having a 
2D system the thickness of the a-Mo2.7Ge film should be 
smaller than ~ 75 nm. 

Going back to the data, it appears that at 100 nm the first 
bending mode becomes relevant, especially for b > 0.15. 
This was not to be expected immediately. For a vortex 
lattice with strong disorder the length at which the first 
bending mode appears might actually be much larger, 
due to the strong renormalization of the elastic constants 
resulting from the forced disorder. Still, the conclusion is 
that apparently this consideration is less important, and 
the scale on which healing starts to appear can be esti- 
mated by considering the first bending mode, while, on 
the other hand, more modes are needed to fully restore 
the VL. 



IV. INTERFACE COUPLING 

A. Experimental data 

In order to investigate how the interface between NbN 
and a-Mo2.7Ge affects the vortex postions, we prepared 
a variety of samples with a NbN layer of 67 nm, a 
50 nm thick a-Mo2.7Ge layer and a ~3.5 nm Au layer, 
where we gradually reduced the interface transparency. 
Next to a reference sample, we first prepared a sample 
for which the NbN layer was exposed to pure O2 gas 
18 mbar for 15 minutes) prior to the deposition of 
the a-Mo2.7Ge/Au layers. The oxidized layer produced 
in this way is very thin; oxidation in air at temperatures 



NOx 1 nmAl203 4 nmAl203 




B = 0.7T; size = 472nm B = 0.7T; size = 473nm B = 0.3T; size = 735nm 



FIG. 7: Autocorrelation maps from vortex lattices mea- 
sured at 4.2 K on samples with different interfaces, 
(left) NbN(67)/NbNO/a-Mo2.7Ge(50)/Au, taken at 0.7 T; 
(middle) NbN(67)/Al2O3(l)/a-Mo2.7Ge(50)/Au, taken at 
0.7 T; (right) NbN(67)/Al2O3(4)/a-Mo2.7Ge(50)/Au, taken 
at 0.3 T. 
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FIG. 8: Comparison of cross- sections from autocor- 
relation maps of vortex lattices measured at 4.2 K 
on samples NbN(67)/Au, NbN(67)/a-Mo2.7Ge(50/Au, 
NbN(67)/NbNO/a-Mo2.7Ge(50)/Au, NbN(67)/Al203(l)/a- 
Mo2.7Ge(50) /Au, and NbN(67) /AI2O3 (4) /a- 

Mo2.7Ge(50)/Au, taken at 0.3 T, at different magnetic 
fields as indicated. In the left-hand panel, the data on 
NbN(67)/Au were taken at 0.2 T, all others at 0.3 T. 



of 180 °C or less is known to oxidize only the first 
one or two atomic layers [8]. To reduce the interface 
transparency even more we deposited an insulating 
AI2O3 interface layer of 1 nm or 4 nm thick in between 
the NbN layer and the a-Mo2.7Ge layer. 
In Figs. EJ [HI we present a selection of the acr- images 
and cross-sections as produced from the raw-data vortex 
images measure at 4.2 K for various values of B. From 
the acr-data it can be clearly seen that on changing 
the interface from 'perfect and thin' to 'insulating and 
thick' the decoupling of the vortex lattices in the two 
superconducting layers rapidly sets in, showing already 
a clear structure for the sample with the oxidized NbN 
surface. From the cross-sections, for the oxided surface 
we find a value for Tc/ao of about 4, for the 1 nm 
AI2O3 about 7, and for the 4 nm AI2O3 larger than 8 
(the maximum for this image size). In other words, the 
4 nm insulating interface has completely decoupled the 
a-Mo2.7Ge vortex lattice from the one in the NbN. 
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B. Shear force versus Josephson coupling 

To gain some qualitative understanding for this be- 
havior, we assume that the vortex positions in the tem- 
plate layer are the result of a competition between two 
forces. One is the Josephson coupling force, which keeps 
the vortices in the template in registry with the those on 
the other side of the interface. The other is the elastic 
force which wants to restore the ordered Abrikosov lat- 
tice, and which we equate to the shear force necessary 
to displace a vortex line from its equilibrium position in 
the lattice. In principle, also the electromagnetic force 
works to maintain registry, but a simple estimate shows 
that the electromagnetic coupling is much smaller than 
the Josephson coupling for the very thin interfaces under 
consideration. We refer this point to the end of the dis- 
cussion. 

Starting with the shear force, for small displacements 
u « ao from the ideal lattice positions (the elastic limit) 
the typical force Fgh necessary to move a single vortex 
line inside the template layer of length dt over a distance 
u is given by 



■ sh 



CQQUdt 



^oBc2(t)udt 
167r/ioA^ 



6(l-0.2%)(l-6)^ (5) 



where the extrapolation formula of Brandt [9] is used for 
the elastic modulus cqq. This approach does not take into 
account that the VL may be disordered (not in equilib- 
rium), and may contain defects which soften the lattice, 
so it should be considered as an upper limit. By rewriting 
C66 = ^§/(6)/(16^/ioA2a2), with/(6) = (l-0.2%)(l-6)2 
and using the fact that h = B/Bc2 ^ ^o/(<^o^c2), we can 
also write Fsh in terms of the quantity eo = ^q/ (47r/ioA^) 
as 



FsH = 



eoudtf{b) 



4al 



(6) 



For the Josephson coupling force, we start from the 
expressions for the Josephson energy Ej as function of a 
displacement of the vortex segment given in ref. [10] 



ttEj 
eodi 



' 0.707 (^)' H^f) , u<2rj 



2.828 



(7) 



1.414 , 2rj <u 



Here, the Josephson length rj is given by rj = jdi, 
with di the thickness of the insulating layer and 7 the 
anisotropy parameter of the system defined by 7 = 
= (mx/^||)^^^ = These expressions for 

Ej are derived from fitting the numerical solutions of 
a model which consist of two stacks of pancake vor- 
tices in an anistropic superconductor (characterised by 
its anisotropy parameter 7), displaced with respect to 
each other by u. For very small displacements, the ex- 
pression for Ej is equivalent to an earlier one given by 
Koshelev and Vinokur : Ej{u) = ^{^)^ H^j/u) 



For larger displacements, Ej becomes linear in ix, which 
means the Josephson force Fj becomes constant and is 
simply given by 



Fj 



2.8 ep 



(8) 



We use this expression as an estimate for the maximum 
Josephson force [12]. Comparing it to the shear force of 
Eq. [6l the condition for Fsh ^ Fj becomes 



u ^ 3.6 ao 

ao ~ 



7 dtf{b)' 



(9) 



Not surprisingly, the result says that the necessary dis- 
placement for the shear force to win from the Josephson 
coupling force is inversely proportional to the length of 
the vortex in the template layer (d^), and also inversely 
proportional to the anisotropy parameter. For our ex- 
periments, typical values are dt = 50 nm and B = 0.3 T, 
which corresponds to ao = 89 nm and f{b) ~ 1, lead- 
ing to {u/ao) ^ 6.4/7. Assuming a value of u/ao ^ 0.3 
for full disorder, the template layer will become ordered 
due to shear stress for 7 > 20. This is a very reasonable 
number for the 4 nm AI2O3 interface layer. In a study on 
multilayers of a-Mo2.7Ge/Ge [13], the following values for 
(di^j) are quoted : (4.5 nm, 7); (5.5 nm, 13), (6.5 nm, 
22). The strong increase with di obviously reflects the 
exponential dependence of the coupling through an in- 
sulating interface, and the numbers given above make 
it quite reasonable to expect that decoupling has not 
yet been achieved for Al2 03(l), but occurs for Al203(4). 
Also, the result indicates that the field dependence of the 
shear modulus is hardly important in comparison with 
the dependence on the interface thickness and for given 
di^ either coupled or decoupled behavior can be expected 
for all fields, as is indeed observed. 

With the expression for the Josephson coupling force, 
it is simple to show that the electromagnetic coupling 
force between two vortex segments can be neglected in 
our considerations. It was shown by Clem [3] that, for 
two superconductors characterized by penetration depths 
Ai,2 = A, with thickness di^2 = d {d <^ X), separated by 
an insulating layer of thickness di {{di <C A), the max- 
imum force as function of the displacement of the 
vortex segments with respect to each other is given by 



87r/io \A^ 



d 



(10) 



Here we neglected a weak field dependence which occurs 
in the regime of interacting vortices, which would lower 
Fem further. Writing Eq. 1101 as -Fem 

eo(i^/(2A^), we 

find that Fem/Fj = -fd^/{2\^). With d of order 50 nm 
and A of order 500 nm for both NbN and a-Mo2.7Ge, Fern 
is an order of magnitude lower than Fj at the thickness di 
where Fsh starts to win from Fj. We should note, how- 
ever, that this may not be the case for superconductors 
with smaller values for A. 
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V. CONCLUSION 

In conclusion, we have used the technique of vor- 
tex core imaging by low-temperature STM to study the 
restoring of initial disorder in a vortex lattice. The disor- 
der was engineered by depositing a weakly pinning super- 
conducting thin film of a-Mo2.7Geon a strongly pinning 
layer of NbN. By varying the thickness of the weakly 
pinning layer, and the interface conditions between both 
layers, we could reach several conclusions with respect to 
the healing mechanisms. For clean interfaces we found 
that the thickness where reordering of the vortex lattice 
is first observed can be linked to the first bending mode 
of the vortex lines, without the need to renormalize the 
elastic constants in the presence of such strong disorder. 
For our a-Mo2.7Ge layers, this thickness is around 75 nm. 



We then used layers with a thickness well below this value 
in order to investigate the restoration of the vortex lat- 
tice when making the interface more insulating. Here we 
find that the results can be understood in a straightfor- 
ward way by comparing the Josephson force working to 
align vortex segments over the interface with the elastic 
restoring forces inside the weakly pinning layer. 
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